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ROTATIONALLY SYMMETRIC SOLUTIONS TO THE CAHN-HILLIARD EQUATION 


ALVARO HERNANDEZ AND MICHAL KOWALCZYK 


Abstract. This paper is devoted to construction of new solutions to the Cahn-Hilliard equation in 
Staring from a Delaunay unduloid Dr with parameter r € (0, r*) we find for each sufficiently small e a 
solution u of this equation which is periodic in the direction of the Xd axis and rotationally symmetric with 
respect to rotations about this axis. The zero level set of u approaches as £ ^ 0 the surface Dr. We use a 
refined version of the Lyapunov-Schmidt reduction method which simplifies very technical aspects of previous 
constructions for similar problems. 


The Cahn-Hilliard equation 


( 1 . 1 ) 


1. Introduction 


Ut =—A{e'^Au — F'{u)) in 

1 ^ = 0 onan, 

ou 

■^ie^Au —F'iu)) = Q on 911, 

ov 


where E is a double-well potential, is a model that describes the process of phase separation of two components 
of a binary alloy. Here 12 C d > 1, is a bounded domain represents the place where the isolation of the 
components takes place, and u, as usual, denotes the outer normal on 912. The function u represents the 
concentration of one of the components and e is the range of intermolecular forces. The double-well potential 
F(u) corresponds to the free energy density at low temperatures, and in this paper we will take 


F{u) = ^ ~ —u. 

From now on we will denote F'(u) = —f{u). 

Equation (1.1) can be derived from the gradient flow of the Helmholtz-free energy functional 


(1.2) E^{u) = J (^F{u{x)) -b ie2|Vu(a;)p^ dx 

in subject to the average concentration to be constant, i.e. 


(1.3) 




where m G [—1,1] (see [16], [17] for details). Note that constant functions u = ±1 are minimizers of this 
functional subject to m = ±1. 
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Stationary solutions of (1.1) satisfy the Euler-Lagrange equation (with f{u) 

e^Au + f{u) = Se in fl, 


-F'{u)) 


(1.4) 


1 

M 


|2 = o 

ov 

/ u dx = m 
Jn 


on do.. 


where 6^ is a Lagrange multiplier. 

Using L-convergence approach Modica [23] showed that minimizers Ue of (1.2) subject to constraint (1.3) 
E-converge, as e ^ 0, to the function 1 — 2xaoj where xAq is the characteristic function of an open set Aq C U. 
Moreover dA^ fl fl is locally a surface of constant mean curvature (CMC surface for short). Geometrically 
the set Aq minimizes the perimeter functional Pern (A) among the sets A C U whose volume is fixed. A 
generalisation of these results was given by Sternberg [26]. Furthermore Hutchinson and Tonegawa [18] studied 
limits of general critical points (1.2) and showed that their limits are locally minimal or CMC surfaces. On 
the other hand it is known [20] that if a set A C H is an isolated mimimizer of the perimeter functional 
subject to the constant volume constraint then there exists a sequence of minimizers itg of ( 1 . 2 ) which F 
converges to A. This result can be used to construct solutions to (1.4) at least in dimension 2 , see [ 8 ]. The 
most complete construction is due to Pacard and Ritore [24] who proved the following: if M is a compact 
Riemannian manifold and A is a non degenerate minimal or CMC sub manifold of M which divides M 
into 2 disjoint components then for all sufficiently small e there exist critical points of ( 1 . 2 ) whose 0 level 
set converges to N. The counterpart of this theory for the time dependent problem (1.1) was developed 
among others by Alikakos, Bates and Chen [1] who proved that as e 0 the time evolution of interfaces is 
governed by the Helle-Shaw problem-of course CMC surfaces are stationary points of the flow. More detailed 
description of the Cahn-Hilliard flow and key spectral tools can be found for instance in [3], [5], [4], [2], [7] 
and the references therein. Additional examples of stationary solutions for the singular perturbation problem 
in a bounded domain have been constructed in [28], [27], [ 6 ]. 

Scaling variables x i-A x/e in (1.4) and letting e tend to 0 leads in a natural way to the following problem: 


(1.5) 


Art +/(u) = (5, in 


In dimension d = 1 there is an obvious solution of this problem when 5 = 0, namely the unique odd and 
monotonically increasing heteroclinic solution H of the ODE, which satisfies 


( 1 . 6 ) 


H" + f{H) = 0, in R, 
i7(±oo) = ± 1 . 


If a € is a unit vector and 6 S R then the function 


u(x) = H(a ■ X + b), 


is also a solution of (1.5) with 5 = 0. When 5 7 ^ 0 there exist radially symmetric solutions to (1.5), see [25]. 
In both cases the level sets of the solutions are CMC surfaces, in the former case their mean curvature is 0 
and in the latter case it is positive number equal to where i?o is the radius of the level set of the 


solution. The radially symmetric solutions in R 
sets are cylinders, which again are CMC surfaces. 

Dilating of the independent variable by a (large) factor e~^ 


can be lifted trivially to R'^ giving solutions whose nodal 


> 0 


e ^x. 


we obtain an equivalent form of (1.5): 

(1.7) eAit H — f(u)=ie, in R . 

e 

where we have denoted | Clearly, if is a solution of (1.7) then v{x) = '«e(ex) is a solution of (1.5). 

On the other hand, if u is a solution of (1.5) then Ue(x) = u(|) is a solution of (1.7). In particular this means 
that while phase transition of the solutions of (1.5) are of order 1, for the solutions of (1.7) they are of order 
£. Thus the latter are more “concentrated”. In the sequel we will focus on solving (1.7). From what we 
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have said above about the singular perturbation problem it is clear that level sets of these solutions should 
converge, as e tends to 0, to CMC surfaces in In fact we expect (on the basis of formal calculations in 
section 2.3) that the Lagrange multiplier 

= f H'{sfds + 0{e), 

^ Jr 

where S is the surface of the phase transition and Hs is its mean curvature. 

We will now introduce a family of embedded CMC surfaces which are good candidates to be the limits of 
the nodal surfaces. We recall that Delaunay unduloids [14], [15] are a one parameter family Dt, t G (0,1) of 
embedded, periodic CMC surfaces of revolution in When the real parameter t tends to 1 the surfaces Dr 
approach the straight cylinder while when r —> 0 they become an array of identical spheres arranged along 
the X'i axis. It turns out that Delaunay surfaces can be constructed in any dimension d > 3 and from now 
on by Dr, t G (0, r*) we will denote the family of Delaunay surfaces in We note that the parameter r, 
satishes: 

2 

^ (d-l)(d- 2 )('^- 2 )- 

Again, in the limit r r* the surfaces Dr approach the straight cylinder, see [19] or Section 2.1 for details. 

It is convenient to ’’normalize” the Delaunay surface and suppose that the mean curvature of Dr is 1 for 
all T G (0,T,). We will also denote by the vector field normal to Dr- Let us notice that the surface Dr 
divides the space into two disjoint components , such that R'^ \ Dr = D+ U 0“, where Nr points towards 
Djl. By changing the orientation of Dr if necessary we can chose Nr in such a way that Djl contains the Xd 
axis. 

Our result is: 

Theorem 1.1. For all r G (0,t,) there exits Er > 0 such that for all e G (0,Et) the problem 

(1.9) eAu H— f(u) = ie in R"^ 

E 

has a solution Ur,e> which is one-periodic in the direction of the Xd-axis and rotationally symmetric with 
respect to rotations about the same axis. As e —>■ 0 we have = 1 + 0{e), and Ur,e satisfies 

Ur^e 1 as e ^ 0 in 
Ur,e —i” — 1 as E —0 in 

uniformly over compacts. 

Remark 1.1. In this paper we took f{u) = u — , which is the standard nonlinearity for the Cahn-Hilliard 

equation. Theorem 1.1 holds for more general nonlinearities of bistable, balanced type, namely f G such 
that f{u) = —F'(u) where F is a double well, even potential with non degenerate wells at ±1. Rather 
straightforward modifications required in the proof of the more general setting are left to the reader. 

We will explain now the implementation of the Lyapunov-Schmidt reduction we used in this paper and 
discuss the differences between our approach and the older implementations which can be found in [24] and 
[11], [12]. Let us first recall the standard Lyapunov-Schmidt reduction method in its abstract version (see [9]). 
Given Banach spaces X, Y and a linear operator A: X ^ Y and a continuous, nonlinear operator N: X ^ Z, 
we are to solve the problem: 

(1.10) Ax-N{x)=0. 

Let 

A/'(A) = Y CX, 7^(A) = W CZ, 

and let iry, 'xw be the projections on the corresponding subspaces. There exists a bounded linear operator 
K ■. W ^ TZ{I — TTy) (the right inverse of A) such that AK = / on IT and KA = I — tty, and moreover the 
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equation ( 1 . 10 ) is equivalent to the equation 

x = y + z, y &Y, z&'R{I — 'Ky) 

(1.11) z - K-KY/Nijj + z) =Q, 

(7 - TTw)N(y + z) = 0. 

In applications the Lyapunov-Schmidt method consists of reducing (1.10) to (1.11), solving the first equation 
for z with y given (which usually can be done by a fixed point argument) and replacing this solution in the 
second equation to obtain the reduced problem 

( 1 . 12 ) {I-TTw)N{y + z{y)) = G. 

In practice several complications may arise and we will illustrate this considering a related to our problem 
which was treated by Pacard and Ritore [24], and in many aspects it is similar to problem we consider in this 
paper. Let A7 be a compact, closed manifold of dimension n and N G M a minimal n — I dimensional sub 
manifold which divides M into two disjoint components. Consider the problem 

(1.13) + u(l — u^) = 0, on M. 

We say that N is non degenerate if the Jacobi operator of N 

Jn = Ajv + + RiCg (vtv, vjv) 

has empty kernel (Ajv is the Laplace-Beltrami operator on N, is the norm of the second fundamental 

form, RiCg is the Ricci tensor on M and vjv is the normal vector to N). The result proven in [24] is: given 
a non degenerate, minimal sub manifold N oi M for each sufficiently small e there exists a solution of 

(1.13) such that the zero level set of itg approaches TV as £ ^ 0. Moreover, converges to ±1 uniformly over 
compacts of the two disjoint components of TV7 \ TV. 

Let us explain now the implementation of the Lyapunov-Schmidt reduction in [24]. It is expected that 
for a; £ TV7 near TV we should have Us{x) = 77(£“^dist (a;, TV)) -|- (p, where dist(-,TV) is the signed distance 
function on T17, H is the unique odd, monotonically increasing solution of —77" = 77(1 — 77^) in R and p 
is a small perturbation. The problem to solve for p amounts to inverting the linearized operator around 
77(£“^dist (a;, TV)) which has form 

L = AM + /'(i7(£-Mist(-,TV))). 

It is known that the norm of L~^ is large due to local translational invariance of the problem. Thus we need 
to perturb TV as well. To describe this perturbation we consider a manifold Nh to be a normal graph over 
TV described by a smooth and small function h: N —R. Furthermoe we let th{x) = dist (a:, Nh) to be the 
signed distance from Nh- Then we look for a solution of the form 



Now both h and p are unknowns. The problem to solve for p is 

Lhp = T{h,p), 

where Lh is the linearized operator around 77^^^. The Lyapunov-Schmidt reduction strategy amounts to 

projection of the above equation onto the function 77'^^^ and its complement, denote this last projection 
by TTh- This leads to a problem for p 

TThLhP = TrhJ'{h,p), 

which we solve first for a given h, and the problem for h 

(1.14) jN,h = gih), 

which we solve next {Jn^, is the Jacobi operator of Nh). Let us discuss (1.14). We notice that the expression of 
in local coordinates will depend in general on h and its derivatives up to order 3, while the Jacobi operator 
is itself only a second order operator. This loss of regularity was dealt with in [24] using a regularisation 
procedure. In a series of papers [11], [12], [13] del Pino, Kowalczyk and Wei introduced a slightly different 
approach to circumvent this problem. It amounts to considering perturbation in the normal direction of the 
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fixed manifold N so that u = + ..., where now t is the signed distance from N and h is a smooth, 

unknown function defined on N. Equation (1.14) takes form 
(1.15) JNh = g{h), 

and the problem of the loss of regularity is thus avoided. The problem is now reduced to finding a fixed 
point of Jjv o G{h), using for example Banach fixed point theorem. To do this we need to know that G is at 
least Lipschitz in h. In both implementations of the Lyapunov-Schidt reduction described above this is rather 
complicated technical point since G depends in a non explicit, non local and non linear way on h. This is 
mainly due to the fact that the linearized operator Lh still depends on h through the potential f 

In this paper we propose still another modification to the method. The idea is simple: instead of working 
with an approximation of the form u = i? +... with h unknown we will improve the initial approxima¬ 

tion to w{t, y) = 77(I) -f ..., t being the signed distance to N and y G N in such a way that we do not need 
to ’’move” N anymore. In other words h will be determined with some sufficient precision before setting up 
the Lyapunov-Schmidt reduction, which with this modification will look like the abstract setting described 
at the beginning. This way we avoid both the loss of regularity and technical difficulties due to complicated 
character of the nonlinear function G^h). This is described in detail in section 3.1. 


2. Preliminaries 


2.1. The surfaces of Delaunay. The Delaunay unduloids Dr, t G (0,t,) are CMC surfaces of revolution 
in Thus for instance in one can parametrize them in the form 

x(t, 0) = {p(t)cos0,p{t)sin0,t), 


where p{t) solves 


Ptt -(1 + Pt) — (1 + — 0- 

P 


However, in this paper we will use mostly isothermal coordinates of Dr C 


( 2 . 1 ) 


Xr{s,e) = i(Te'^^(^)0,K,(s)), (s,e) eKx 


where functions {ur, Hr) are the unique solutions of the following system of ODEs: 


( 2 . 2 ) 


(5«n,)2 + ^r2(e--+e(2-dW)2=o, a.a,(0) = 0, n,(0) < 0, 

dsKr - ^^2(6'^- + ^ 0, R^(0) = 0. 


We will now summarize some basic facts about the Delaunay surfaces and their isothermal parametrization. 
The function CTt is periodic, and consequently the surfaces Dr are one-periodic along the x^-axis: namely if 
2Tr denotes the minimal period then 

Dr = Dr -f 2Tr&d- 


Clearly we have the relation 

Tr = ^Kr{2Sr), 

where 2sr is the minimal period of Ur- 

The Jacobi operator Jr of Dr is defined by: 

(2.3) Jr:=/\Dr+\Ar\^, 

where is the Laplace-Beltrami operator on Dr and \Ar\'^ is the square of the norm of the second 

fundamental form of Dr- The Jacobi operator is of fundamental importance in this paper and to understand 
well its properties we will first consider the special case 7 = 3. In the isothermal coordinates (s,0) € M x 5^ 
its expression is given by: 

(2.4) Lr = + dg + COSh(2cr^) } . 

The Jacobi fields on Dr, which are elements of the kernel of Lr are of three types: 
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(1) Jacobi fields arising from infinitesimal translations. For any e G |e| = 1 the constant Killing 
field associated to translations 

X I —> e 

induces the following Jacobi fields 


= e-Nr, 

where N^. is the unit normal vector to 0^. The coordinate vectors e^, j = 1,2,3 generate three 
linearly independent Jacobi fields corresponding to translations of Dt in the directions of the 

coordinate axis. We note that in the isothermal coordinates 


$T.e 3 ^ ^ j = 1 , 2 . 

T e - 

It is important to notice that the Jacobi fields $r’ ^ are bounded. 

(2) Jacobi fields arising from infinitesimal rotations. Let e G |e| = 1 be such that e • 63 = 0. The 
Killing vector field corresponding to the rotation about the vector e is: 

X I—S’ (x • e)e3 - (x ■ 63)6. 

We define the Jacobi field associated to this vector field by: 

= [(x • 0)03 - (x • 03)0] • Nr. 

There are clearly two linearly independent Jacobi fields associated to the rotations. They are: 

and they correspond to rotations about the coordinate axis. Note that in isothermal coordinates 
functions , j = 1 , 2 grow linearly as functions of s. 

(3) Jacobi field associated with the variation of the Delaunay parameter. We define: 

= -drXr ■ Nr. 


This Jacobi field is somewhat harder to write explicitly however it can be shown that the function 
<i>^(s) is linearly growing. 

In summary, the Jacobi operator Lr has at least 6 explicit Jacobi fields which are either linearly growing or 
bounded. By a result of Mazzeo and Pacard [22] we know that these are all Jacobi fields with at most linear 
growth. To explain this let us observe that by separation of variables the equation Jrtp = 0 separates into a 
sequence of problems 

Lrjip = 0, Lr,j = {d'^ + cosh{2ar) - f} , |j|=0,1,.... 

Then we have: 


Proposition 2.1 ([22]). The homogeneous problem Lrjif = 0 has the following solutions: 

(1) one periodic and one linearly growing solution when j = 0 or \j\ = 1; 

(2) two solutions (pfj{s) which satisfy: 

+ Sr) = 

with 

^r.j Ct,j ^ 0 ; 

when jjj > 2. The numbers fr,j “re the indicial roots of the operators Lrj and correspond to the 
behavior of the solutions of the homogeneous problem at ±00. 

This basic facts can be generalized for the Jacobi operator of Delaunay surfaces in d > 3. We will 
summarize them now and refer the reader to [19] for details. We have the following at most linearly growing 
Jacobi fields: 

(1) The are d bounded, periodic Jacobi fields arising from infinitesimal translations. They will be denoted 
by j = l,...,d. 
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(2) There are d — I Jacobi fields arising from infinitesimal rotations in the direction of xj axis, j = 
1,..., d — 1. We will denote them by ^, j = 1,..., d — 1. and they correspond to rotations about 
the coordinate axis. Note that in isothermal coordinates functions , grow linearly as functions 
of s. 

(3) There is one Jacobi field associated with the variation of the Delaunay parameter <I)^ = —drX-r ■ Nr, 
and it is linearly growing. 

2.2. Fermi coordinates and shifted Fermi coordinates near a CMC surface. Let S be an embedded 
CMC surface in and let denote its mean curvature. By N we will denote its unit normal. We will 
assume that there exists a tubular neighborhood Afs of S of width 25 in which we can introduce local system 
of coordinates (Fermi coordinates) {y,z) G S x {—5,5) by setting: 

X I— {y, z), where x = y + zN{y). 

We suppose that this map, which we denote by Y, is a diffeomorphism from Ms to T' x {—5,5) whenever 5 is 
taken sufficiently small. In the sequel we will use the inverse of this map 

Y-^-. Sx {-5,5)^ Ms 

{y,z) I—^ X. 

Given a function w. Ms we define its pullback Y*w to 17 x {—5, 5) by the diffeomorphism Y as: 

Y*w{y,z) =woY~^{y,z). 

For technical reasons we will chose later the size of the tubular neighbourhood 5 depending on e but for now 
on we just take 5 small. 

Next we will define shifted Fermi coordinates. To do this we let /i: 17 —>■ R be a given smooth function 
such that the map 

XI —> {y,t), where x = y+{t + h{y))N{y), 

is a diffeomorphism from Ms into 27 x (—<5, 5). We will denote this map by Yh and by we will denote its 
inverse, finally by Y^w we will denote the pullback of w: Ms —>■ R"^ by Yi,: 

Yhw{y,t) = woY^^{y,t). 

It will be convenient to have at hand expressions for the Laplacian in Fermi and shifted Fermi coordinates. 
To derive them by 27^ we will denote the surface 27 + zN i.e. the original surface 27 shifted in the direction 
of the normal by z. Locally near 27 we have 

A = 

We denote zM^^z = A^;^ — A^;. The operator Ms,z is a second order differential operator. To expand the 
curvature term we use the well known formula: 

= = Hs + z\As? + + 0{z^) = Hs + z\As\^ + z'^Qs,.- 

j=l -5 j = l 

where kj are the principal curvatures of 27. In summary we have: 

A = As+d'^ — {Hs + z\As\^)dz + zE>s,z + 

The reason we expanded the Laplacian in this way will become clear later on. From this it is easy to obtain 
a formula for the Laplacian in the shifted Fermi coordinates: 

(2.1) A = Ax: + (1 + |V£'/ip)9^ — {Hx + Axh + {t + h)\Ax\^)dt + (t + h)Mx,t+h + (t + 7i)^Qi;,t+/i- 

Anticipating the content of the next section we introduce the stretched shifted Fermi coordinate 

t 
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Formal consideration will show that an approximate solution We of the Cahn-Hilliard can be obtained if we 
assume that it is a function of the form: 

yhWe{y,t) = H +o(l), 
where H is the heteroclinic solution of (1.6) 

As before we have a diffeomorphism and its inverse '■ S x (—|, j-) Ms, and for any function 
w: A /5 —?► we define its pullback by by: 

Taking onto account formula (2.1) we get 

(2.2) A = Aj; + e ^(1 +1Vi;/i|^)9.^ — e ^(i?i; + Ax'h+(£t + h)|A5;|^)9t + (£t; + /i)B£^et+/j + (£t + /i)^Qx’,Et+/!,. 


2.3. Formal expansion of the solution of the Cahn-Hilliard equation concentrating on E. For the 

purpose of formal calculations we will assume that the the solution of (1.7) near A is a function ui, which 
depends on the stretched and shifted Fermi coordinates (y,t), in the following way 

(2.3) = t/(t)+£Vo(y,t),, 

for some functions U and ipQ which we will determine. Moreover, we will assume that 

h = e^ho, 


where hp is a constant to be chosen. 

To determine U and tpo we write the error N^(w) — '■= eAw + ^w(l — w^) — ig in the local coordinates 

N,{w) - 4 = s{e-^d^U - e-^Hsd^U + e-'^f{U) - £-14 

+ dtV'o - eHEdttjjo + f'(U)'>po - (t + eho)\As\'^dtU 

(2.4) +lf(^u;)--f{U)-ef{U)% 

e e 

+ s'^As'tpO — £^(t + £/lo)|A£'|^clt'0o) 

+ ((£t + h)lBs^gt+h + (st + h)‘^Qx:,et+h) (U + £^'^0)}. 

In order to get as small as possible this approximation we have to get rid of the first three terms of the right 
hand side of the expression above, since they show the lower powers in e. To write things compactly let: 

Soiw) := d'^w - eHsdtW + f{w), 

Lqw := d^w — eHsdtW + f'{U)w. 

With this notation and the ansatz (2.3) we can write the problem in the form: 

So{U) + e'^ LQ^l)Q + Qo{U + s’^ipo) = s^E- 

where Q{U + e'^ipo) represents the rest of the terms in (2.4), and we also have to determine the Lagrange 
multiplier 

Thus we take the Lagrange multiplier to be a number such that the following ODE has a unique, 
monotonically increasing solution U 

So{U)=U"- eHsU'+ f{U)=e£e, in R, 

fiU{±^)) = £4. 

This function is easy to find by perturbing the heteroclinic solution H. 

Note that since / is odd symmetric we have f7(±oo) = ±1 + where 

/(±1 + CTe ) = £ 4 . 

Also, we have 

i^=iQ-\-0{e), 4 = f H'{s)‘^ds, 

^ Jr 
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and 

Uit) = H{t) + 0{e). 

Next, we will determine the 0{e^) correction ipQ. Ignoring terms of order 0{e^) we get the following 
equation to solve: 

(2.6) -eHsdtilie + (1 - 3C/^)V'o = {t + eho)\As\'^dtU. 

It is convenient to consider, more generally, an ODE (with the right hand side possibly depending on the 
variable y) of the form: 

(2.7) dl^p - eHsdtip + (1 - 3U^)ip = g(y, t). 

A solution of this problem can be found by the variation of parameters formula. Indeed, the fundamental set 
of the ODE is spanned by the functions 

(2.8) StD = O((cosht)’'*), tE(t) = C)((cosht)‘'^), t ^ ±oo. 

with 

(2.9) ^ - Y^-4t(±oo) + ^ (eHs + Y^-4t(±oo) + , 

and 

i(oo) = 1 — 3(1 + de)^, i( —oo) = 1 — 3( —1 + fJe)^. 

We can assume that the Wronskian at 0 is 1. If the right hand side of (2.7) satisfies 


[ 9 {'^,y)dtU{t)e = Vy e A 

Jr 


( 2 . 10 ) 

we can write ip = Q{g) where 

(2.11) Q{v){t,,y) =—dtU{t) f W{s)e~^^^^v{s,y) ds + W{t) ( dgU{s)e~^^^’^v{s,y) ds. 

Jo J — oo 

Note that the orthogonality condition (2.10) guarantees that the function Q{g) is exponentially decaying 
whenever g is exponentially decaying (in t). To be more precise let us assume for instance that 

| 5 (y,t)|(cosht)^ < C, 

with /i £ (ry + eHs, —ly], where g = maxjry"'', 77 “} <0. Then we have 

|V3(y,t)|(cosht)'" < C, 

as well. 

Using this we can determine the function i/jq. To do this we need to chose ho such that (2.10) is satisfied, 
in other words: 

hence 

With this choice we define: 

Note that we have: 
and as a consequence 

IV'o(t,y)|(cosht)^ < C, 0 < fi <-g. 

Sometimes it is convenient to derive a more refined estimate taking into account the fact that the leading order 
term on the right hand side is tdtU{t) = 0(|t|(cosht)’* )), |t| ±00. Thus, we consider (2.7) assuming 

that 

( 2 . 12 ) 


f {t + eho)dtU{t)^e dt = 0, 
Jr 


hn = 


= o{i), 


° e dt 

'0o(y,t) = g{{t+eho)dtU)\As\'^. 

(|t| + e|/io|)|9tl7(t)|(cosht)'" < C, 0 < g < -g, 


l5(y,t)|(l + |t|)^(cosht)'" < C, 
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with ^ G (ji + sHs, —r]], and /3 S K. By a simple argument we have as well: 

(2.13) |C/(g)(t,y)|(l + |t|)^(cosht)'^<C. 

3. Delaunay solutions of the Cahn-Hilliard equation 

3.1. The Lyapunov-Schmidt reduction. While our formal considerations in the above section were valid 
for any embedded CMC surface S in in what follows we will focus on a special example when S = Dt, 
i.e. it is a Delaunay unduloid. Since we are interested in functions which are periodic in the direction of the 
Xd axis with the minimal period equal to that of Dt we will introduce the manifold Dt which is obtained by 
identifying the set Dt C {xd = 0} with the set Dt Cl {xd = 2Tt}. The set Dt is homeomorphic to the d — 1 
dimensional torus 

Fist we note that the approximate solution w = U + e'^ipo is so far only defined in Afg, which is a tubular 
neighborhood of Dt- To extend w to the whole space let us define 


IHI(x) = 



if X e l)+, 

if X e Dt , 


where D^, Dt denote, respectively, the interior and the exterior of Dt- 

Let us notice that the function w approaches HI exponentially. Indeed, we have 

|i;>(y,t)-F4H(y,t)| y G Dt, ^€(^-^,0, 

for any /i € (0, — ry). To make the definition of w* precise we let y to be a cutoff function such that x(s) = 1, 
when |s| < ^ £Hid x(s) = 0, when |s| > 1. Next, we define a cutoff function y* supported in Ms by: 

Ye*hX*i^) = X (y) • 

We can define the global approximate solution w* by 

(3.1) w*{x) = u;(x)y*(x) + H(x) [l - y*(x)]. 

Now we look for a solution of the equation (1.7) in the form 

u = w* + (f 

where (y? is a small (in a way to be specified) function. Thus our problem can be stated: find tp: 

R, which is one-periodic with period 2 Tt in the Xd direction, such that 

(3.2) Ng{w* + ip) = ie, in 
where S^Tt is a circle of radius 2 Tt, 


^ ^ X S2Xt 


jd-l „ cfl 
5. X 02Tt , 


Ne{u) = eAu -G ^14(1 — M) 


and if, is the Lagrange multiplier defined in (2.5). Let us recall that we want our solution to be rotationally 
symmetric. That is, if by TT-g we denote the rotation of R'^ about the Xd axis by angle d then we should have: 

{w* + ip)(px) = (w* + ip){7l0x). 

Since we already have (by definition) 

w*{x) = w*(TZ 0 x), 

than as a result we will have ip{x) = ip{TZ 0 x) as well, as can be seen easily from the proceeding construction. 

Since the function ip appearing in (3.2) is expected to be small it is natural to expand the nonlinear operator 
IVe and write: 

L^ip = -Ns{w*) - Qei^p) + 4 

where 

Leif = DNg{w*)(p, Qei^p) = N{w* + (p) - Ng{w*) - Lgip. 

The strategy, based on the Lyapunov-Schmidt reduction is clear. Indeed, we expect that due to the d 
dimensional bounded (and periodic) kernel of the Jacobi operator JV which is associated to translations of 
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Dr in the directions of the coordinate axis ej, j = 1,..., d, the linear operator Lg should have d dimensional 

T e 

kernel spanned, roughly speaking, by the functions Zr’e ^ where 

(3-3) Y*,,Z^p {y,t) j = l,...,d, V{y,t) = dtw = dtU{t) + e'^dtipo{y,t). 

Notice also that in general any function Z^ such that 

Ye*h^e{y,t) = $(y)V(y,t), 

is also “almost” in the kernel of in the sense that Y*f^{L^Zg) = o(l). We introduce a linear subspace of 
L^{Dr X R) of functions that are orthogonal to Z^ by: 

(3.4) X = e L'^{Dr xR) 

By n we denote the orthogonal projection on X. We set ip = where 

Yl^pW e d:', Yl^p^ = zvex^. 

Finally we split our problem into two equations: 

(3.5) n o Y*f,[N,iw* + + P^) - 4] = 0, 

(3.6) (Id - n) o Ylf^[N,{w* + (^11 + P^) - 4] = 0. 

When solving (3.5) we use the fact that the associated linear operator is coercive on X. To solve (3.6) we 
will make use of the theory of solvability of the Jacobi operator Jf^ . An additional, somewhat technical, 
step which we have omitted in this informal discussion is to ’’transfer” the original problem from the space 
of functions defined on R'^"^ x S 2 T to a space of functions defined on Dr x R. We will explain these details 
in section 3.3 but first we will introduce and study a linear operator which is essentially the expression of 
in the Fermi coordinates of Dr- 



3.2. Linear theory for a model problem. In this section we will develop the necessary theory to deal 
with the operator L^. To this end we will consider the operator 

(3.7) he -(^Dr +£(t+ £^o)x(et/J)|A^j2)at+e“V'('*^), 

where x(s) is a cutoff function supported in (—1,1) and equal to I in (—1/2,1/2). Note that this operator 
is defined for functions </: Dr x R —>• R (and not just functions defined on Dr x (—-, -)). It is clear that 
Y*^Le ~ Le. Although the function w = t/ + e^^o depends on both variables (y, t) in some sense the operator 
he separates variables. To see this, with dtW = dt{U + e'^ipo) = V, we consider functions of the form: 

p{y,t) = V(y,t)Z(y). 

Observe that, by construction oi w = U + e'^ipo, combining equations (2.5) and (2.6) multiplied by we get 
e~^d^w - + e{t + ehQ)\Ap,^\'^)dtw + e~^f{w) = 4 - e^(t + £ho)\Ap,J^)dtipo + e~^Q{e'ipo), 

where Q{v) = f{U + w) — f{U) — f'{U)v. Differentiating this equation in t we get for V = dtW. 
e + e(t + eho)\Af)^\‘^)dt'V - e\Ap,j'^V + e ^/'(w)V =|^9t4o 

- e^(t + eho)\AfjJ‘^d^-ipo - e~'^dtQie‘^ipo)- 

From this, using the definition of Lg in (3.7) we get: 

(3.9) L4:^) = he(VZ) = SVjf,^Z + Be{Z) 

where is the Jacobi operator on Dr and 


Be{Z) =+ eho){l - x*)\Ap I^ZdtV- 


“h eZ V — E 


2iE^^ V • ^ f) Z 


(t + eho)\Afj^ l^d^ipo + \ Aj^^ 


Dr I 




£-1 [f{w)drw - f{U)drW - e^f"{U)d,U^o] Z. 
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We note that 


(3.10) i?,(^) = O(e')ll^llci(D0- 

Identity (3.8) and its consequence (3.9) is the key calculation which allows to use the usual Lyapunov-Schmidt 
reduction scheme, as w explained in the introduction. Indeed, if we had taken as the approximate solution 
only the function U then differentiating the equation (2.5) for U' = dtU we would have gotten 

e-^d^U' - + e-^f'{U)U' = 0 . 


This equation, unlike (3.8), does not carry any information about the geometry of Dr besides its mean 
curvature which is constant. Following the method of [24] or [11], [12], [13] we would have to perturb the 
surface Dr additionally introducing new unknown functions in our problem. With the approach presented 
here this is no longer necessary and the Lyapunov-Schmidt procedure in this version is simpler. Recalling 
that the linearization of the mean curvature operator is the Jacobi operator which depends on the second 
fundamental form, we see that the operator is naturally compatible with the geometric context of our 
problem. To put it differently: the operator Lg is, up to negligible terms, the correct linearization of the 
Cahn-Hilliard operator near a solution whose zero level set is the constant curvature surface Dr- 

To develop invertibility theory for Lg we will we employ two basic facts. First, we observe that on the 
subspace: 

y '■= t) = V(y, t)Z(y) I Z{Y)<^>r’^^ (y) dy = 0, j = 1,..., , 

we have 


(Lg(/j,(/?) > 


Second, when we consider ip G X (space X is defined in (3.4)) and g G L^{Dr x R) such that is a bounded 
solution of the problem 


L,p = g, Lg=eA^^-he ^d'^^+e 


then we have 


(3.11) 11^1! 

L^{DrXR) < Ce\\g\\ L^{DrXR)- 

To prove this estimate we use a contradiction argument which relies on the fact that from 

[ Kp-/'(id)u2 >Cllul|i.(R), if / vm'it)dt = 0, 

Jr Jr 

it follows that the bilinear form 

Bg((^) := {Lep,(p), 

is coercive on X. In the same way as (3.II) is shown one can prove: 

(3.12) ll‘f’lliLi(DTXR) — ^II5 |Il2(^^xR)- 

The loss of the factor of e in this last estimate is related with the scaling of the problem. We refer the reader 
to [11] or [13] where results similar to estimates (3.11) and (3.12) were proven. 

At the same time we can use (3.12) and the coercivity of the bilinear form Bg((/j) to solve the equation 

(3.13) nA'Lg(/j = g, 


where 11;^ is the projection on A, in X. To do this we write Lg = Lg -|- (Lg — Lg) and use a perturbation 
argument. The solution will still satisfy estimate (3.12). The perturbation argument is as follows: for g S A 
we solve 

Let/) = 5 + cH', 

where c = . Then we define a map 


Gx{g) =(j)-V 


f(/>V dt 
JWdt' 
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Note that Gx: X ^ X. Next we check 


Wnxh^Gxig) - g\\ 

l2(£i^xR) <o(l)||5ll 

Therefore Hxl'eGx is invertible as a map from X to itself and we can define 

=Gx{l^eGx)-\ 

Moreover it is rather straightforward to show that a solution to (3.13) will satisfy estimate (3.12). 

We will use these observations to solve the following model equation: 

(3.14) Le(p = 5(y,t), 

where we will assume initially that g £ LP{Dr x R). We look for a solution in the form ip = + (p-*-, where 

pW g A", = zv e x^. 


We write 

= na’Le((/5ll + (p-^) + + V^"^), 

and then we need to solve 

nA-LeVsIl = Wxg - Uxl^eg}^ 

= ^X^g - 

The idea is that terms and are of smaller order because Lgl/ = o(l) so that the coupling 

between the two equations is rather weak. Another important point is that 

=eVJtj^Z + B,{Z), 

where B^{Z) is small (see (3.10)). We decompose accordingly g = g-^ ^ g-^ = SV and Bs{Z) = 

Be{Z) + B^{Z), B^{Z) = Te(Z)V and look for a solution of the system: 

n^L,(^ll=gll-i?|(Z), 

(3 15) f 

eJo Z + T,(Z) = S- / (L,^ll)Vdt + Vc,$pT 

Note that in the second equation we have introduced Lagrange multipliers Cj to be determined. 

Although the two equations in (3.15) are coupled but this coupling is weak and we can solve the system 
without any difficulty using invertibility of Tlx and jTp , a fixed point argument and estimates (3.10), (3.11) 
and (3.12). We leave this to the reader. 

Given that we can solve (3.15) our purpose is to find suitable estimates for the solution of the problem 

= 5(y,t) 


on X assuming that 

5 (y,t) = C>(e"'"l’^l), |t|-:>oo. 

In particular we would like to know that :/9(t,y) = C>(e“^l’^l) as well. This is straightforward by comparison 
principle once we know for example that p is bounded. Thus the main issue is to obtain L°° control for p. 
We will go a little further now and show how to control a priori certain weighted Holder norms of p^^ and p'^ 
(see decomposition in (3.15)). 

In order to simplify the argument and avoid keeping track of negative powers of e appearing on the right 
hand side of various estimates we will rescale the y variable. Thus we will introduce: 


y = 


y 

e 


t = t. 
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We will denote Dr,e = ^-Dr- We consider the manifold Dr^e x R equipped with the product metric and the 
associated Levi-Civita connection. We dehne weighted Holder norms on DT,e x M: 

ll“llc“'“(D.eXR) = SUp(cOsht)^||u||^o..(i^^x(t-lA+l))> 

t£R 

ll^llc/.’“(D^,eXR) = ll“llc”’“(DT,eXR) + II^D^,eXR“llc”’“(D^,eXR)’ 

ll^llc 2 ’“(D^.eXR) = ll“llc”’“(DT,eXR) + II^DT,eXR“llc”’“(D^.eXR) + II^D^,eXR“llc“’“(D^.eXR)- 
Above and respectively denote the gradient and the Hessian on the manifold Dr,e x M. 

Given functions u and g on : Dr x R and we introduce functions 

uCy.'t) = w{Y,i) =’wiey^i), 5(y,t) = £g(ey,t). 


and set 


Left = + d~u + f'{w)u + eqdiu, 


where 

9 := + e(t + eho)x{ei/S)\Af^^ (ey)p 

is a bounded function. The linear problem we consider is: 

(3.17) Axl^eU = g, in Dr,e x R, 

where now we assume simply that the right hand side satisfies the orthogonality condition (c.f. (2.11)): 


/ 5(y,'t)V(ey,t)dt = 0. 
Jr 


For this problem we can derive a priori estimates using the method of [11], [12], [13]. From this, by a 
perturbation argument, we will be able get corresponding estimates for satisfying the first equation in 
(3.15). Let us explain briefly the main steps. 

Step 1. We consider a problem of the form: 


(3.18) Ay(l) + d^(l} += inR'^-^xR. 

The following result is known ([11], [12]): 

Lemma 3.1. Let (f) be a bounded solution of (3.18). Then cj) = cH'(t), with some constant c. 


Step 2. Consider now equation (3.17) and assume that g € C^’°'(Dt,£ x R), with p S (0, jjyj). The we have: 

Lemma 3.2. There exists a constant C > 0 such that for all sufficiently small e any bounded solution of 
(3.17) satisfies: 

ll^llc2’“(D^,eXR) — ^ll5llc°’“(D^,eXR)- 

A proof of this lemma, which relies on Lemma 3.1 and a contradiction argument, follows the same lines as 
the proof of Lemma 5.2 in [13] (see also similar results in [11], [12]). 

Now we should go back to the original variables. We define weighted Holder norms on Dr x R similarly as 
in (3.16): 




(3.20) 


“(D^xl 




We note that if for a given function u: Dr x 


Dr 


Dr 


i''^llc“’“(D^x][ 

{^llc“’“(D^x][ 


I^D,xr'“IIc[;’“(D,x][ 


(3.21) 




e-\\d'(D^uU^^ + 


we set M(y, t) = u(£y,t) then we have 


^ Dr^''C'^(rirXM.) - 

0<k-\-m<l 0</c+m<^ 




Dr >a,tJ,,Dr xR’ 
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where weighted Holder seminorm. Consequently by £f/°‘{DT x 

functions on Dr x R with the norm 

With this definition we have 

while with the notation of (3.21) we get 


we denote the space of 


c- 




From this and Lemma 3.2 it follows for ^ = 0,1, 2: 
(3.22) 


“ll£^’“(T)xXR) - Il5ll£);'“(ij^xl[ 


norm of in (3.15) 


The procedure described above shows that we can control the size of £j^°‘{Dr x 
obtaining: 

(3-23) llf^’“(i3^xR) — (IIS’''I|£2 ’“(£i^xR) + ^ llff ll£2’“(i’xXR))- 

Finally, we notice that for the second equation in (3.15) using elliptic theory we can get Holder estimates 
and since (p"*" = ZV we find: 

(3-24) llv? IIc^’“(_d^xR) — II 5 llc“'“(i!)^xR)' 

3.3. The linear problem in the whole space. Now we will use the theory outlined above to solve the 
following problem: 


(3.25) 


£^^+-f'{w*)ip = g{yi), in 
£ 


pd—1 


X S, 


2Tr ) 


From what we have said above it is in general not possible to find a solution with a reasonably bounded norm 
unless the right hand side satisfies some extra conditions, or equivalently, we need to introduce Lagrange 
multipliers that correspond to natural invariances of the problem. Thus, we will solve 

1 

_ 4-1 ( 

£‘ 


(3.26) 

where 


£l£ip + -f'{w*)ip = g{-x.) + x*^^]Zr;e'^ inR"' ^ x. S2Tr 

i=i 


^e%X*(t) = x(et/(5), y;)^Z^^®^(y,t) = V(y,t)$p^(y), j = l,...,d 

The idea is to solve (3.26) by gluing a solution defined near Dr and another one defined away from Dr- To 
describe this construction rigorously we need some preparation. We introduce the function g(x) as follows: 


g(x) = 


/'(I + CTe), dist (x,!)^) > 5/2, 
/'(-I + CTe), dist (x, I?^) <-5/2, 


and otherwise q(x) is a smooth function such that min{/'(l + + CTe)} < 9(x) < max{/'(l + 

+ <^e)}- Note that q{x) = —2 + 0(e). Finally, we need another cutoff function x such that 
XX* = X* (take for instance F/*^x(t) = x(et/25) and chose 5 smaller so that the Fermi coordinates are 
defined in A//^). We want to find a solution of (3.26) in the form ip = x*'P ° Ye,h + i/, where the function p 
solves: 

d 

(3.27) heP = (x{9 + X*'^'^jZr;e' + “ L^)ip - [x* , L^]ip + £~'^[q - f {w*)]'tp}y in Dr X R, 


J=1 


and the function 1 / solves 


(3.28) eAiP + e ^[(l - x*)f{w*) + x*q]^ = (I - x *){9 + X*^CjZ^p - [x*,L^]ip}, 

i=i 


in ^ X S2Tt ■ 
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It is clear that multiplying (3.27) by x* and adding the equations (3.27)-(3.28) and using the fact that 
XX* = X* gsl the solution to our problem. In the above and in what follows we abuse slightly notation 
writing for instance as a function defined on Dr xM and as a function defined on x S 2 Tr ■ It is understood 
that in the latter case we take o To avoid complicated notions we will omit the composition with 1)./j 
or Y~^ whenever it does not cause confusion. Thus the commutator [x*,L^\(p in x S 2 Tr is 

[x*,Ls]if> = 2eVipoY^^h^x* + SV’°Ye,h^X*, 

while in Dr x R we have to first express in local coordinate (y,t) (written as Y*i^Lg) and calculate 

[x*,Y*HLe]^. 

The function g on the right hand side of this equation satisfies the following general assumptions on its 
asymptotic behaviour: 

Widx II£ 0 ’“(d^xr) - 

(3.29) II(5X*)^IIc«-(d.xR)<C’ 

In addition we assume that g is rotationally symmetric about the Xd axis, namely if by TZg "we denote the 
rotation of about the Xd axis by angle 9 then g{TZgx) = g{x). 

In order to solve this coupled system we need to make sure that all terms on the right hand side that 
involve tp and are small in suitable weighted Holder and Holder norms respectively. It is at this point that 
we need to chose the parameter S in the definition of the tubular neighbourhood Afg small and dependent on 
e. Thus we take (5(e) = e^/^. This means in particular that for x G Ms we have et(x) = 0(e^/^). For reasons 
that will become clear soon we will also chose the Holder exponent a in the definition of C^’‘^{Dr x R) and 
C'^’°‘{Dr X R) to be in the interval (0, j^). Finally, the parameter p will be always taken in the interval (0, Iryj). 
Considering equation (3.28) we have the following: 


Lemma 3.3. If u is a solution of 

1 , 


eAu - [(l-X*)/'(w*) + x*g]u = 5 , in 

e 


bd—1 


X S2Tr 


then we have an a priori estimate: 

(3.30) l|u||cc<.(R<i-ixS2T^) ^ 




||co.“(H‘i~ixS2T^)' 


The proof of this lemma is straightforward and it is omitted, for similar results see for instance [10]. From 
this we get readily an a priori estimate for (3.28): 

d 

(3.31) lll/'||cc<»(Rrf-ixS2T^) ^ ^ “11(1 - X*)ffllc‘>'“(R‘^~ixS2T^) + 6 (T^ 1*^11 + ll5P|lch°r_D.,xR'i)- 

i=i 

From the theory developed in the previous section we can also obtain an a priori estimate for (3.27). If we 
write 

d 

0 = x{5 + X*X! + (^^ “ “ [x*Me]p + e~^[q - f'{w*)]ili}, 


then we have using (3.23) 
(3.32) 


£m’“(DtXR) — 


< Cei-“ 


^“'“(D^xR) 


+ e^ll0^llc«" 


(D^xl 


and using (3.24): 

(3.33) II^^IIc^.“(d.xR) ^ ^£"'llfl^llc2'“(D.xR)- 

We note that the weighted norms we use for and are scaled differently with e. This slight nuisance is 
a result of our choice of the original scaling of the Cahn-Hilliard equation. We observe as well that with our 

definitions jj • ||fO.“(£)^xR) ~ II ' IIc°-“(DtxR)' 
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We will now estimate ||0''||£o.“(i3^xR)- observe that, with ^ = ^11+ we have: 


+ £ 


||(x(Le Le)(f) ||c°'“(i)^xR) - " ll‘^" Ilf2’“(-b^ xR) 

ll(x[X i^e]‘P) ll£“’“(_D^xR) — ll£,i'“(Z>^xR) + 

||(x£”M9-/'(w^*)]V')"Lo,.(£,^^r) < C'e-dl^llc«.“(R-i-ixS2T.)- 
Next we estimate the orthogonal complement of these functions 

-■^-)‘^)^L«'“(T).xR) ^ ^(II<^"||£2'“(D.xR) 


^ llc2’“(b,xR) 

-ce-l/^ll,x-L 


,) 


)||(^ ||gl,<»(^^xR))i 


IIc2'“(£>xX1 


l) 




ll(x[x ll£“’“(i3^xR) — )(ll‘^" ll£L“(i3^xR) ll*^ llci’“(£lTXR) 

||(xe“^[9- < C'e“dl' 0 llco.“(R‘i-ixS 2 Tj 

We can estimate the parameters cj by projection of g onto Zr,’e^ ■ Using the above estimate we get: 

|Cil < uj IKx g) + II (x(L£ - Le):^) ||g0,<»(j5^x"" 


II (x[x% ^e]^) ^ + II (x£ ^ [q - f'iw*)]i’) ^ Ilc«-(i3^ XR) } 

< C” |||(x 5 )‘'‘IIc“’“(Z)tXR) + ll‘^dl£2’“(bxXR) + ^^lli^'''llc2’“(li>^xR) 

~^0{e )(ll‘^dl£L“(i;)^xR) 11''^ llcJ’“(i!)^xR)) ^ II'0 ||cO'“(R‘^-1xS2t^)| 

^IIV'IIc°>“(R‘*-1xS2t^)} 


< c 


- £ IIc2’“(£>xXR) ^ 


ll£i’“(i:)^xR) + I 

'{iKx 5)‘'‘llc;;’“(b^xR) + II‘^"II£ 2 ’“(-d^xi[ 

Now we use estimates (3.32)-(3.33). After rearranging terms suitably and using e^““(5“^(e) = o(l) to absorb 
(^11 in the first inequality below we get 

Ile 2 ’“(-DTXR) - {IKx ff)''llf 2 '“(L)xXR) IK^ Ilc;i’“(i3^> 

(3.34) 


- e dll/'llco.“(R'^-ixS 2 T,) + '^(^)II'^'''IIc2’“(£'xxr)}’ 

' " ' " "f2’“(i3^xR) + ^ lll/'llco.° 


V)}- 


ll'i^ IIc 2 -“(£)^xR) - (IKx 5) llc;i’“(£)^xR) + ll‘^" llf 2 ’“(_D^xR) + ^ IIUIICO.“(R‘i-i XS2T, 

From (3.31) we get as well for £ = 0,1, 2: 

(3.35) 

IIV'llcL°(Rrf-ixS2T^) ^ ^ “{||(1 - X*)ffllc'>'“(R3) 

+ e (IIXffllc“’“(b^xR) + ll‘^''li£^'“(b^xR) + Il‘^''"llc2'“(f)^xll 

Using the fact that (5(£:)e““ = o(l) to absorb term <5(£)||:^-*-||g2,c^^^^gj appearing on the right hand side of 
the hrst inequality in (3.34) and combining these estimates we get 

(3.36) 


■ lll/’llcO.°(R-i-lxS2T^))}- 


£^’“(ji)^xR) - 


IIc2'“(L)xXI1 


< 


Ce^ “ {ll(X5)"llf:«.“(b^xR)+£ ^'5(e)||(x £/)-^llc“’“(b 
+£ ^ “^(e)ll(l “ X*)ffllco.“(R‘i-ixS2T^)} > 
Ce-^ {ll(X'7)^llc«.“(i3.xR) + e^”“l(x <?) 


l-au~ „1|||| 


ll£2’“(.bxXl 


-2-a 


11(1 - X*)ffllc“.“(R'i-lxS2T^)} > 


and 

(3.37) 


IIV'II 


- 1/8 


C2> 


-lxS 2 T,) ^ C'e ^ “11(1 - x*)ff||co.“(R‘i-ixS 2 T^) + ^*(6' ■ JIIXl/llc2'“(ii,xi 


, , , ~ ^ , .. —/ 

Using these a priori estimates we can solve the system (3.27)-(3.28) by a standard fixed point argument. To 
do this we replace the functions (/H, (p'^,'4’ on the right hand side of the system by known functions tF 

which satisfy estimates of the same type as (3.36)-(3.37) but with constants bigger that those appearing in 
(3.36)-(3.37). Then we have a map 
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from a certain ball in the space x K) x C^’°‘{Dt x K) x x S 2 T^) into itself. This and the 

Lipschitz character of this map being evident from the way we have derived a priori estimates allows for an 
application of the Banach fixed point theorem. We leave the details to the reader and simply state this as a 
result for (3.26). 

Lemma 3.4. For each sujficiently small e there exists a solution of (3.26) in the form ip = x*^ ° ^s.h + V' 
such that estimates (3.36)-(3.37) hold. 


3.4. Proof of Theorem 1.1. Now we can finish solving the nonlinear problem 

(3.38) L,p = i,-N,{w*)-Q,ip). 

As we saw above we need to modify this equation by introducing Lagrange multipliers. Thus we will consider: 

d 

(3.39) L,p = 4 - N,{w*) - Qe{p) + X*Y. 

i=i 

To solve this problem we use a fixed point argument and the linear theory in Lemma 3.4 above. The first 
task is to calculate the size of the error of the approximation 4 ~ This is straightforward using the 

definition of w* and formula (2.2). We recall here that h = e^ho, where ho is a constant and consequently 
this last formula simplifies significantly. We can write: 


4 - N,{W*) = x *[4 - N,{w)] + [N{w*) - x*N{w) - (1 - x*)N{M)] = + A 2 , 

since 4 = A^£(IH) in supp (1 — y*). Using exponential decay oi w — (±1 + a^) when t -A ±00 we get easily: 

11(1 - X*)^2||co.“(R'i-ixS2T,) ^ CQe~^^ . 

To estimate Ai some standard calculations which we will omit are needed (c.f Section 2.3). As a result we 
get 

(3-41) IlK.hX^ll] llc)l'“(D.r xR) — , 

11(1 - X*)^i|IcO’“(R‘*-1xS2t^) ^ C'oe , 

where Cq, and 9 are positive constants. Now we use the linear theory developed in the previous section to 
solve the nonlinear problem (3.39). Thus we write p = x*F ° + '4’y further decompose p = pi^ -\- p-^ 

where S A n C‘fp°‘{br x R), p-^ &y C\ Cfp°‘{br x R) and ip G C^’“(R^“^ x S2Tb- To set up a fixed point 
scheme we fix functions p ^^, p-^ and ip in these sets such that 

llf:2’“(A^xR) — ’ 

(3-42) \\p llcA“(D.rxR) 4 ATe , 

ll^llcA<»(R<i-ixS2T^) 4 ATe”®® , 

where AT is a large constant to be chosen and 9 G (9/2,9) is a constant. Let us denote the right hand side 
of (3.39) by g. It is evident that under the assumptions (3.42), and with a suitable choice of the constants 
a > 0 and p G (0, |? 7 |) we can solve the problem (3.39) for functions (p'^/p-^,if) which again satisfy (3.42). 
Thus we have a non-linear map 

(F'^F^y^) '-1 (p'/p^,1p), 

of this set into itself. To show that this map is a contraction is straightforward, using the quadratic nature 
of the nonlinear function Q(p). At the end we have a solution of the problem: 
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where Zr,’e^ is the (approximate) element of the kernel of the linear operator associated with translation 
in the direction of the Xj axis, see (3.3). To show that in fact 

Cj = 0, j = 

we need: 


Lemma 3.5 (Balancing formula). Let X = o,jdxj he the infinitesimal generator of translations or rotations 
in For any C^(M'^) function it holds: 


(3.44) 


div 


£ 1.-^ i2 1 


-IVmM - -F{u) X(u) - eXluWu 
2 e 


= — [eAu+ -F'{u)]X{u). 


We will take Xj = for some 1 < j < d and integrate the balancing formula over the cylinder = 
Br X S2Tx- Using (3.43) and Green’s theorem we get: 


[ (^\Vu\^ —-F{u) + Uj dS — f dx^udnUdS 

JaCn V2 e J JacR 


'Cr 


X*^3’ZT’e'' I dxjU. 

^t' = l 


Ir 


Hr 


IIIr 


The first integral Ir is 0 on the top and the bottom of Cr and on the other hand, using the asymptotic 
behavior of the solution we get finally lim/j_>oo d/j = 0. In the second integral the integrals over the top and 
the bottom of Cr cancel because u is periodic. Then, from exponential decay of the derivatives of u we get 
lim/{_>.oo Hr = 0. Finally, we note that 

hdxj u 

hence 


r,€ ’ 


IIIr = c 


/. IZr’e'' f + o(l) Cj', 

JDxXS. j,^^ 


from which we get immediately c^- = 0, j = 1,..., d. 


References 

[1] N. D. Alikakos, P. W. Bates, and X. Chen. Convergence of the Cahn-Hilliard equation to the Hele-Shaw model. Arch. 
Rational Mech. Anal, 128(2):165-205, 1994. 

[2] N. D. Alikakos, X. Chen, and G. Fusco. Motion of a droplet by surface tension along the boundary. Calc. Var. Partial 
Differential Equations, ll(3):233-305, 2000. 

[3] N. D. Alikakos and G. Fusco. The spectrum of the Cahn-Hilliard operator for generic interface in higher space dimensions. 
Indiana Univ. Math. J., 42(2):637-674, 1993. 

[4] N. D. Alikakos and G. Fusco. Slow dynamics for the Cahn-Hilliard equation in higher space dimensions: the motion of 
bubbles. Arch. Rational Mech. Anal., 141(1):1-61, 1998. 

[5] N. D. Alikakos, G. Fusco, and V. Stefanopoulos. Gritical spectrum and stability of interfaces for a class of reaction-diffusion 
equations. J. Differential Equations, 126(1):106-167, 1996. 

[6] P. W. Bates and G. Fusco. Equilibria with many nuclei for the Gahn-Hilliard equation. J. Differential Equations, 160(2):283- 
356, 2000. 

[7] X. Chen. Spectrum for the Allen-Cahn, Cahn-Hilliard, and phase-field equations for generic interfaces. Comm. Partial 
Differential Equations, 19(7-8): 1371-1395, 1994. 

[8] X. Chen and M. Kowalczyk. Existence of equilibria for the Cahn-Hilliard equation via local minimizers of the perimeter. 
Comm. Partial Differential Equations, 21(7-8):1207-1233, 1996. 

[9] S. N. Chow and J. K. Hale. Methods of bifurcation theory, volume 251 of Crundlehren der Mathematischen Wissenschaften 
[Fundamental Principles of Mathematical Science]. Springer-Verlag, New York-Berlin, 1982. 

[10] M. del Pino, M. Kowalczyk, F. Pacard, and J. Wei. The Toda system and multiple-end solutions of autonomous planar 
elliptic problems. Adv. Math., 224(4):1462-1516, 2010. 

[11] M. del Pino, M. Kowalczyk, and J. Wei. On De Giorgi’s conjecture in dimension TV > 9. Ann. of Math. (2), 174(3):1485-1569, 
2011. 

[12] M. del Pino, M. Kowalczyk, and J. Wei. Entire solutions of the Allen-Gahn equation and complete embedded minimal 
surfaces of finite total curvature in Journ. Diff. Ceometry, 93(1):67-131, 2013. 

[13] M. del Pino, M. Kowalczyk, and J. Wei. Traveling waves with multiple and non-convex fronts for a bistable semilinear 
parabolic equation. Comm. Pure Appl. Math., LXVI:481-547, 2013. 

[14] G. Delaunay. Sur la surface de revolution dont la courbure moyenna est constante. J. Math. Pures Appl., 6:309—320, 1841. 

[15] J. Eells. The surfaces of Delaunay. Math. Intelligencer, 9(l):53-57, 1987. 







20 


ALVARO HERNANDEZ AND MICHAL KOWALCZYK 


[16] P. C. Fife. Models for phase separation and their mathematics. Electron. J. Differential Equations., pages No. 48, 26 pp. 
(electronic), 2000. 

[17] P. C. Fife. Pattern formation in gradient systems. In Handbook of dynamical systems, Vol. 2, pages 677-722. North-Holland, 
Amsterdam, 2002. 

[18] J. E. Hutchinson and Y. Tonegawa. Convergence of phase interfaces in the van der Waals-Cahn-Hilliard theory. Calc. Var. 
Partial Differential Equations, 10(1):49—84, 2000. 

[19] M. Jleli. End-to-end gluing of constant mean curvature hypersurfaces. Ann. Fac. Sci. Toulouse Math. (6), 18(4):717-737, 
2009. 

[20] R. V. Kohn and P. Sternberg. Local minimisers and singular perturbations. Proc. Roy. Soc. Edinburgh Sect. A, lll(l-2):69- 
84, 1989. 

[21] M. Kowalczyk, Y. Liu, and F. Pacard. The space of four ended solutions to the Allen-Cahn equation on the plane. Ann. 
Inst. H. Poincare Anal. Non Lineaire, 29(5):761-781, 2012. 

[22] R. Mazzeo and F. Pacard. Bifurcating nodoids. In Topology and geometry: commemorating SISTAG, volume 314 of Con- 
temp. Math., pages 169-186. Amer. Math. Soc., Providence, RI, 2002. 

[23] L. Modica. The gradient theory of phase transitions and the minimal interface criterion. Arch. Rational Mech. Anal., 
98(2):123-142, 1987. 

[24] F. Pacard and M. Ritore. From constant mean curvature hypersurfaces to the gradient theory of phase transitions. J. 
Differential Geom., 64(3):359—423, 2003. 

[25] L. A. Peletier and J. Serrin. Uniqueness of positive solutions of semilinear equations in R’^. Arch. Rational Mech. Anal., 
81(2):181-197, 1983. 

[26] P. Sternberg. The effect of a singular perturbation on nonconvex variational problems. Arch. Rational Mech. Anal., 
101(3):209-260, 1988. 

[27] J. Wei and M. Winter. On the stationary Cahn-Hilliard equation: bubble solutions. SIAM J. Math. Anal., 29(6):1492-1518 
(electronic), 1998. 

[28] J. Wei and M. Winter. Stationary solutions for the Cahn-Hilliard equation. Ann. Inst. H. Poincare Anal. Non Lineaire, 
15(4):459-492, 1998. 

Departamento de Ingenieria Matematica and Centro de Modelamiento Matematico (UMI 2807 CNRS), Universi- 

DAD DE Chile, Casilla 170 Correo 3, Santiago, Chile. 

E-mail address: cdiernandezOdiin.uchile.cl 

Departamento de Ingenieria Matematica and Centro de Modelamiento Matematico (UMI 2807 CNRS), Universi- 

DAD DE Chile, Casilla 170 Correo 3, Santiago, Chile. 

E-mail address: kowalczyOdim.uchile.cl 


